Survival of a diffusing particle in an expanding cage 
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We consider a Brownian particle, with diffusion constant D, moving inside an expanding d- 
dimensional sphere whose surface is an absorbing boundary for the particle. The sphere has initial 
radius Lo and expands at a constant rate c. We calculate the joint probability density, p(r, t\ro), 
that the particle survives until time t, and is at a distance r from the centre of the sphere, given that 
it started at a distance ro from the centre. The asymptotic (t — > oo) probability, Q, obtained by 
integrating over all final positions, that the particle survives, starting from the centre of the sphere, 
is given by Q = [4/I> + l)A" +1 ] b n exp[-(a^) 2 /A], where A = cL /D, b n = (c£) 2v /{J v+1 (a?)] 2 , 
v = id — 2)/2, and a™ is the nth positive zero of the Bessel function J v (z). The cases d = 1 and 
d = 3 are especially simple, and may be solved elegantly using backward Fokker-Planck methods. 



S 1 

o 

£ 

i 

-t— > 

•4-3 



o 
o 



o 
in 
q 

in 
o 

o 



X 



I. INTRODUCTION 

First-passage problems for stochastic systems have at- 
tracted renewed interest in recent years [1], notably in 
systems with infinitely many coupled degrees of freedom, 
where the effective stochastic process describing a single 
degree of freedom is often non-Mar kovian However, 
even in systems with few degrees of freedom, and Marko- 
vian dynamics, problems with moving boundary condi- 
tions have proved difficult to solve. 

In this paper we consider a single diffusing particle, or 
Brownian walker, moving within a d-dimensional sphere, 
of initial radius Lo, which is expanding at a constant rate 
c. The surface of the sphere defines an absorbing bound- 
ary for the particle. We present exact results, for any 
space dimension d, for the probability that the particle 
survives up to time t. 

In an earlier paper Q this problem was solved for the 
special case d — 1, in the limit t — *■ oo, using an elegant 
method based on the backward Fokker-Planck equation. 
In section [II] we recap how this method works and show 
how it can be extended to the case d = 3. The method 
does not, however, appear to be useful for other values 
of d. Why this should be so becomes clear later in the 
paper, where we obtain a complete solution for general 
d. 

To our knowledge, the results presented here (and in 
3) are the first exact results for this type of problem. 
Approximate methods have been developed by Krapivsky 
and Redner [4J in the limit of slow ( "adiabatic approxi- 
mation" ) and fast ( "fast approximation" ) motion of the 
absorbing boundary. The relevant dimensionless param- 
eter is A = cLq/D, where D is the diffusion constant of 
the particle. The slow and fast limits correspond to small 
and large A respectively. 

In section IIIII we show that the adiabatic approxima- 
tion of Krapivsky and Redner can be modified to obtain 
an exact solution of the usual (forward) Fokker-Planck 
equation for the survival probability of the particle for 
any time t, and any initial position within the sphere, for 
any space dimension d. 

The outline of the paper is as follows. In section HT1 
we present the backward Fokker-Planck equation for the 



infinite-time survival probability, and its solutions for 
d = 1 and d = 3. In section HTT1 we show how the forward 
Fokker-Planck equation can be solved to obtain the sur- 
vival probability in any space dimension, for any time t 
and for an arbitrary starting point within the sphere. In 
section HVl we extract the infinite-time survival probabil- 
ity. The previous results for d = 1 and d = 3 are recov- 
ered as special cases. For general d, it is not straight- 
forward to extract analytically, from the exact solution, 
the limiting behaviour for large A, so in section[V]we em- 
ploy the "fast approximation" of Krapivsky and Redner 
to obtain an approximate solution that becomes exact in 
this limit. Section IVll concludes with a brief summary of 
the results. 



II. INFINITE TIME SURVIVAL PROBABILITY 
OF A PARTICLE IN A EXPANDING CAGE: THE 
BACKWARD FOKKER-PLANCK METHOD 

A. Solution in One Dimension 

In we considered a diffusing particle obeying the 
Langevin equation X(t) — r)(t), X(0) — x, where 
77(f) is Gaussian white noise with mean zero and time 
correlator (j](ti)j](t2)) = 2D5(t\ — t 2 ), bounded by a 
linearly-expanding absorbing cage with edges located at 
±(Lq + ct). After a time At the boundaries have moved 
to positions ±(Lo + c(t + At), and the particle has moved 
to x + Ax, where (Ax) = 0, and ({Ax) 2 ) = 2D At. The 
probability, Q(x, Lo,t) that the particle still survives af- 
ter a time t satisfies the obvious equation Q{x, L^,t) = 
(Q(x + Ax, Lq + cAt, t — At)). Expanding to first order 
in At yields the backward Fokker-Planck equation 



dQ _ D d 2 Q c dQ 
dt dx 2 8Lq 



(1) 



The novelty in this approach resides in treating Lq, which 
gives the initial positions of the boundaries, as an addi- 
tional independent variable. 

Discarding the time-derivative term, to directly ob- 
tain the infinite-time limit, and introducing dimension- 
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less variables y = cx/D and A = cLq/D gives the simpli- 
fied equation 



d 2 Q | dQ 
dy 2 dX 



0, 



(2) 



subject to the absorbing boundary conditions 
Q(±A,A) = 0, and with Q(y, A -> oo) = 1. The 
solution is 

oo 

Q(y,X) = J2 (-l)"cosh(ny)e-" 2A , (3) 



n— — oo 



which satisfies both the differential equation and the 
boundary conditions. 

For a particle starting at the origin (y — 0), the sur- 
vival probability Q(0, A) is given by 

oo 

Q(0,A) = E (-l) n e-" 2A ~f-2e-\ A -» oo. (4) 



n— — oo 



We can extract the small-A behaviour by using the 
Poisson sum formula, 



E /( n )= E f( 2 ^ 



n— — oo 



where /(fc) is the Fourier transform of the the function 
/(n), to recast Eq. ([3]) in the form 



— oo 

7 



(5) 



k= — oo 



Solution in Three Dimensions 



We now extend this result to the case of a diffusing par- 
ticle in three spatial dimensions bounded by an linearly- 
expanding, absorbing sphere of radius L(t) = Lq + ct 
using the same backward Fokker-Planck method. 

For general spatial dimensionality d, equation ([2]) has 
the obvious generalisation 



v 2 g + f = 

OA 



(6) 



where the dimensionless spatial coordinate is r = cr /D, 
and A = cLq/D as before, where ro is the initial location 
of the particle within the sphere. Exploiting the spatial 
isotropy, we infer that Q depends on r only through its 
magnitude, r — |r|, giving, for d = 3, 



d 2 Q 2 dQ dQ 



0, 



dr 2 r dr d\ 
This equation has separable solutions of the form 



(7) 



characterised by an index k, from which a general solu- 
tion can be constructed by superposition. Note that solu- 
tions of the form [cosh(Ax)/r] exp(— fc 2 A) are rejected as 
they are not regular at the origin. At this point we make 
no assumptions about the values of k. We find, however, 
that no simple superposition of solutions Qk(r, X) satis- 
fies the required boundary condition Q(X, X) = 0. We 
note, however, that any derivative of Qfc(r, A) with re- 
spect to k is also a solution of Eq. ([7|) since it is simply a 
superposition of two values of k which are infinitesimally 
close together. We therefore try taking the first deriva- 
tive of this solution and postulate a sum of the functions 
dQk{r,X) / dk with integer values of k, in analogy with 
the one-dimensional case, to obtain 



Q(r,A)= E 



cosh(nr) sinh(rtr) 

r 



If we choose the amplitudes to be a,k — 1 for all k, we see 
by inspection that the boundary conditions Q(A, A) = 
and Q(r, A — > oo) = 1 are satisfied. The survival proba- 
bility is therefore given by 



Q(r, A) = — [rcosh(nr) — 2nAsinh(r 



(8) 

For a particle starting at the origin we take the limit 
r — > 0, giving 



Q(0,A) = E (l-2Afc 2 ) e" 



k=-c 



1-2(2A- l)c~ A , A 



(9) 
(10) 



To obtain an expression suitable for extracting the be- 
haviour at small A, we take the Poisson transform of this 
sum, giving 



2tt 5 / 2 °° 

Q^) = %jt E 

k— — oc 



k 2 e~ 



2tt 5 / 2 

A3/ 2 ' 



(11) 



for A -► 0. 

For general space dimension d, the backward Fokker- 
Planck method does not seem to be useful. The reason 
for this will become clear in the following section. 



III. SURVIVAL PROBABILITY IN A 
d— DIMENSIONAL EXPANDING SPHERE 

We approach the problem in general dimension by solv- 
ing the forward Fokker-Planck equation for the prob- 
ability density, p(r,i|r o ,0), defined as the probability 
density that the particle, starting at position ro within 
the sphere of radius Lq, still survives (has not yet 
reached the absorbing boundary) at time t, and is cur- 
rently at position r. It satisfies the partial differential 
equation dp/dt = D\7 2 p. The boundary condition is 
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p(r, i|ro,0) = for |r| — Lo + ct, and the initial condition 
isp(r,0|r ,0)=<5 rf (r-r ). 

In generalised polar coordinates, the Fokker-Planck 
equation reads 



dt V dr 2 



d 2 p d — 1 dp 



dr 



£-nP 



(12) 



where £n is a generalised angular derivative operator. 
From the spherical symmetry of the problem, p(r, i|i"o, 0) 
only depends on r = |r|, tq = \tq\, and the angle between 
r and i"o. The problem can be simplified by choosing the 
direction of ro as the principal polar axis, and integrating 
out the angular degrees of freedom by defining 



p(r, t) 



1 



dfip(r,t|r O) 0), 



where dQ is an element of solid angle, S d = 2n d / 2 /T(d/2) 
is the integral over the solid angle, i.e. Sd is the surface 
area of the unit sphere in d dimensions, and the depen- 
dence of p(r, t) on ro is implicit. 

The differential equation for p reads 



The initial condition is 



d 2 p d — 1 dp 



di 



dr 



p(r,0) 



5(i 



(13) 



(14) 



Our method of solution for the moving boundary prob- 
lem is motivated by the solution for a fixed absorbing 
boundary at r = Lq. Separable solutions of Eq. (fl"3]) . 
regular at the origin, have the form 



Pk(r, t) 



Jyjkr) k 2 Dt 
e , 



where J u (x) is the Bessel function of order v, and v — 
(d — 2) /2. The absorbing boundary condition, p(L , t) = 
0, selects a discrete set of k- values, k n = a™/L , where 
a" is the nth positive zero of J u {x), to give the discrete 
set of solutions 



P 



n ( r , t ) = M<r/L ) c - (a?i) i Dt/L i 



(15) 



Our trial solution of Eq. (TIB")) replaces L by the time- 
dependent L(t) = Lo + ct, and t/L 2 by f* dt'/L 2 (t'). We 
also multiply the static solution by an unknown function 
of r and t, to give 



P{r,t) = g(r,t)- 



exp 



dt' 



( 16 ) 

At first sight it seems that we have simply replaced one 
unknown function, p(r,t), by another, g(r,t). This is, of 
course, true. However, the resulting equation for g(r,t) 



simplifies greatly for the case of interest, namely L{t) = 
Lq + ct, when it can be exactly solved. The same trial 
solution, Eq. (|16[) . with n — 1 (the lowest mode) and 
g(r,t) — * g{t), was used by Krapivsky and Redner [i[ as 
a general technique for obtaining an approximate solution 
in the adiabatic limit where the wall is moving slowly. 

On substituting Eq. into Eq. (fl~3"|) we obtain the 
following equation for g(r,t): 



-0, (17) 



where dots indicate time derivatives and J' u {x) = dJ v /dx. 
We seek a solution in which the terms that involve Bessel 
functions, and the terms that do not, separately vanish. 
Such a solution exists if the two equations 




g = D (d rr g + -d r g 



d r g = - 



Lr 
2DL 



(18) 
(19) 



are both satisfied. Eq. (fT9"]) can be integrated immedi- 
ately to give 



g(r,t) = A{t) exp 



Lr 2 
4DL 



(20) 



where A(t) is an arbitrary function. Substituting this 
result into Eq. (TT8"|) gives 



.4 



Lr 2 
4DL 



A. 



(21) 



But A(t) depends only on t, and not on r. For a con- 
sistent solution, therefore, we require L = 0, i.e. the 
absorbing boundary must move at constant speed. Solv- 
ing Eq. (|2ip for this case gives 



A{t) = K/L(t) 



(22) 



where K is an arbitrary constant. 

The general solution for the case L(t) = Lo + ct is 
obtained as an arbitrary superposition of separable solu- 
tions, combining equations (fTB|) . (I20[) and (f2"2"| : 



p(r,t) = 



J v 



L(t)r» \L(t) 



L Q L(t) ADL{t) 



(23) 



the summation being over all positive zeros of J v (x) . 

As usual, the amplitudes a n are determined by the 
initial condition, Eq. (fl4|) . exploiting the orthogonality 
property of Bessel functions: 



rdr J,, 



Lo 



J, 



T 2 
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The exact solution is thus 



p(r,t) = 



J, 



- I a™r \ j I a"r \ 
" \L(t)J J »{ L ) 



L L(t)(rr o yS d ^ [-WO] 2 



L L{t) 4D ITTtT L 



(24) 



Eq. (f2"4"]l represents our most complete result. Note 
that p(r, t) depends on the constants ro, c and D 
as well as r and t. To make contact with our earlier 
results for d = 1 and d = 3, we now compute the infinite- 
time survival probability, for given starting radius ro, by 
letting t — > oo and integrating over the final coordinate 



IV. THE INFINITE-TIME SURVIVAL 
PROBABILITY 

To find the probability that the particle survives for 
infinite time, we integrate the probability density over the 
whole domain. We introduce the dimensionless variables, 
p — cro/D and A = cLq/D, in terms of which the infinite- 
time survival probability can be expressed: 

2e P 2 /4A ~jj!%£) 

Q(p, A) = lim 6 , . V „ \ \ 

a '-In{t) (25) 



xc 



where 



Mt) r d-i d f a n r\ -£rg 

= / — ^(iS))^' (26) 



To extract the large-i behaviour, we change variables to 
x = r/L(t) in the integral. Then for large t the inte- 
gral is dominated by small values of x and the Bessel 
function can be replaced by its small- argument form 
J v (z) ~ (z/2) u /T{v + 1), to give, for t ->■ oo, 



/n \d/2 

I n (t) ~ 2^ (f ) 



L(t) . (27) 
Putting this result into Eq. (|25| gives our final result 



Q{p ' X) -yx e h [J, +1 «W e ' (28) 

where we recall that p — cro/D, A = cLq/D, and ^ = 
d/2 - 1. 

It is interesting to consider the special case of a particle 
starting at the origin. Taking the limit p — > in Eq. 
gives 



Q(0, A) 



n\2v 



K) 



r(i/ + l)A"+i^ 1 [J H - 1 (a!J)] 



(29) 



The limiting form for small-A is given by the first term 
in the sum. 

We now compare our results with those obtained using 
backward Fokker-Planck methods in section [TTJ The case 
d = 1 and d = 3 correspond to v = —1/2 and v = 1/2 
respectively. The functions J-i/2(.z) and J\/2{z) have 
positive zeros at z = (2n— 1)tt/2 and z = nir respectively, 
while J 1 2 /2 [(2n - 1)tt/2] = 4/[?r 2 (2n - 1)] and J 2 /2 [n7r] = 

2/(7r 2 n). It is the fact that the Bessel functions zeros are 
uniformly spaced in d = 1 and <i = 3 that makes these 
cases especially simple. 

For d = 1, the final result, Eq. (|29|) . can be written as 



Q(p, A) 



I — OO 

\/^ Yl cos((2n-l)7rp/2A) 



-(2n-l) 2 7r 2 /4A 



n— — oo 

2^2 



(30) 



while for d = 3 it takes the form 



A) 



- P P 2 /4A 



xe 



— -= nsin(n7rp/A) 
PvA ft= __ 



(31) 



For the special case p = 0, these results reduce to Eqs.© 
and pip . One can show this correspondence holds for 
general p by using the Poisson summation formula (with 
p — > y in d = 1 and p — > r in <i = 3) to transform Eq. (|28j) 
into Eqs. ([3]) and ([5]) for e? = 1 and d = 3 respectively. 

It is not simple to use the Poisson summation formula 
on Eq. (|28|) for general d, since we do not have explicit 
expressions for a™ except for v = ±1/2. This means 
that it is not straightforward to analytically extract the 
large A limit of Q(0, A). In the final part of this paper, 
therefore, we extend the "fast approximation" of ref. [H 
to general dimension d in order to investigate the large-A 
behaviour. 



V. APPROXIMATE SOLUTION FOR FAST 
EXPANSION (LARGE A) 

For a rapidly expanding cage in d dimensions we fol- 
low the approximate method used in 0]. This approach 
consists of equating the loss of survival probability to the 
outward probability flux at the boundary. For a particle 
starting at the origin, we define Q(t) to be its survival 
probability at time t. If the cage is expanding rapidly 
the probability flux through the boundary is small and 
we can use the free form for the diffusion propagator: 



p(r,t) = QQt) 



ex.p(-r 2 /4Dt) 
{^Dt) d / 2 



where r is the distance from the origin at time t. We 
equate the rate of decrease of Q to the flux through the 
boundary: 



dQ_ 

dt 



S d L{t) d -"Dd rP {r,t)\ r=L{t) , 
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VI. SUMMARY 



lnQ(i) 



r(d/2)(4D)«*/a 7o 

(L + cu) 2 



x exp 



4£>it 



(32) 



Setting t = oo, we can evaluate the integral for large 
A using the method of steepest descents to obtain the 
infinite-time result 

nt \ ( 2 V^ A ^ -a^ 
Q(po) — exp | e 1 



r(d/2) 

^ , — . d— i 

_ 20FA^ —a 

r(d/2) 



A — > oo, (33) 



which, for the cases d = 1 and d = 3, agrees with the 
large-A limits given in Eqs Igj) and (fTTJ)) respectively. 



In this paper we have derived exact results, in arbi- 
trary space dimension d, for the survival probability of a 
particle diffusing inside a uniformly expanding cage that 
acts an an absorbing boundary for the particle. We found 
very simple forms for the survival probability in one and 
three dimensions by using a backward Fokker-Planck ap- 
proach. Our solution for general d, however, using for- 
ward Fokker-Planck methods, indicates that the simplic- 
ity of the solutions in one and three dimensions is fortu- 
itous. Our method of solution seems to be restricted to 
the case where the boundary moves at constant speed, 
though it will be interesting to pursue the question of 
whether there are any other soluble cases. 
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